# Model (1) in the paper 
#
# Specify the model, setting effects separately for each dependent network
# network effects
# density (in- and out-degree)
multiPtaEff <- setEffect(multiPtaEff, density, name="posPtas", fix=TRUE)
multiPtaEff <- setEffect(multiPtaEff, density, name="negPtas", fix=TRUE)
multiPtaEff <- includeEffects(multiPtaEff, inPopSqrt, name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, inPopSqrt, name="negPtas")
# triadic closure
multiPtaEff <- setEffect(multiPtaEff, transTriads, name="posPtas", include=FALSE)
multiPtaEff <- setEffect(multiPtaEff, transTriads, name="negPtas", include=FALSE)
# balance
multiPtaEff <- includeEffects(multiPtaEff, balance, name="posPtas", include=TRUE)
multiPtaEff <- includeEffects(multiPtaEff, balance, name="negPtas", include=TRUE)
#
# covariate effects
# constant dyadic covariates
multiPtaEff <- includeEffects(multiPtaEff, X, interaction1="Distance", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, X, interaction1="Distance", name="negPtas")
# time-varying dyadic covariates
multiPtaEff <- includeEffects(multiPtaEff, X, interaction1="Trade", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, X, interaction1="Trade", name="negPtas")
#
# income variables
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="CapIncome", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, simX, interaction1="CapIncome", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="Income", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, simX, interaction1="Income", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="CapIncome", name="negPtas")
multiPtaEff <- includeEffects(multiPtaEff, simX, interaction1="CapIncome", name="negPtas")
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="Income", name="negPtas")
multiPtaEff <- includeEffects(multiPtaEff, simX, interaction1="Income", name="negPtas")
#
# trade in services as share of GDP
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="servicesTrade", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="servicesTrade", name="negPtas")
#
# 
# stock of US FDI in the services sector in the country in the year
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="usFdi", name="posPtas")
multiPtaEff <- includeEffects(multiPtaEff, altX, interaction1="usFdi", name="negPtas")
